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ABSTRACT 

We report the discovery of a radio pulsar with a characteristic age of 7300 years, making it one of the 10 
apparently youngest Galactic pulsars known. PSR J 1357-6429, with a spin period of P = 166 ms and spin-down 
luminosity of 3.1 x lO^^ergss"^ was detected during the Parkes multibeam survey of the Galactic plane. We 
have measured a large rotational glitch in this pulsar, with AP/P = -2.4 x 10"^, similar in magnitude to those 
experienced occasionally by the Vela pulsar At a nominal distance of only ~ 2.5 kpc, based on the measured 
free electron column density of 127 cm"^ pc and the electron distribution model of Cordes & Lazio, this may 
be, after the Crab, the nearest very young pulsar known. The pulsar is located near the radio supernova remnant 
candidate G309.8-2.6. 
Subject headings: ISM: individual (G309.8-2.6) — pulsars: individual (PSR J1357-6429) 
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1. INTRODUCTION 

A supernova (SN) explodes in the Galaxy every '^ 100 years 
(see, e.g., Cappellaro 2003, 2004, and references therein). As 
a result, young neutron stars in their various guises are rare, 
as are young pulsars. Nevertheless, the payoff resulting from 
their study can be large, in areas as varied as the physics of 
core collapse, the internal structure of neutron stars, and mag- 
netospheric emission processes. Also, the youngest neutron 
stars are often embedded in compact nebulae powered by rel- 
ativistic pulsar winds or otherwise interact with their host su- 
pernova remnants — as such they can make for magnificent 
probes of their immediate environment. Young pulsars are 
also particularly useful for establishing reliable birth rates of 
this important branch of outcomes of core-collapse SNe. For 
these reasons, substantial effort continues to be devoted to the 
detection of the youngest neutron stars, and three-quarters of 
the Galactic pulsars known with an age r < lOkyr have been 
discovered in the past seven years at radio and X-ray wave- 
lengths. 

In this Letter we announce the discovery of PSR J 1357- 
6429, a very young, relatively energetic and nearby pulsar, 
present its rotational history during a span of 4.5 years that 
includes a large glitch, and comment on the immediate envi- 
ronment of the pulsar and in particular on a nearby supernova 
remnant candidate. 

2. DISCOVERY AND OBSERVATIONS OF PSR J1357-6429 

PSR J1357-6429, with a period P = 166 ms, was discov- 
ered on 1999 October 7 in data collected during the course 
of the Parkes multibeam survey of the Galactic plane (e.g., 
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[Manchester et alJ 1200 lY This survey employs a 13-beam 
low-noise receiver system at a central radio frequency of 
1374 MHz with a recorded bandwidth of 288 MHz. The large 
area covered by one pointing enables the 35 min-long indi- 
vidual integrations that, together with the high instantaneous 
sensitivity, result in the good limiting flux density of about 
0.2 mJy for long period pulsars over a large area along the 
Galactic plane (\b\ < 5°;260° < I < 50°). I n turn this has led 
to the discovery of more than 700 pulsars ([Ma nchester et al 
[200!| iMorris et alJ [2002; Kramer eTaD l2003t iHobbs'etar 
^004^. 

After discovery, as with every newly detected pulsar, we 
began regular timing observations of PSR J 1357-6429 at 
Parkes, using the filterbank-based observing system em- 
ployed in the survey. Typically this consists of the record- 
ing to magnetic tape of raw data for about 15 min in order to 
obtain off-line a pulse profile from which we derive the time- 
of-arrival (TOA) of a fiducial point on the profile by cross- 
correlation with a high signal-to-noise ratio template. The 
average pulse shape consists of a single approximately sym- 
metric peak with FWHM = 15 ms (0.09P). In this manner we 
obtained 125 TOAs spanning the MJD range 51458-53104, a 
period of4.5yr. 

Using the initial set of TOAs together with the TEMPO 
timing software*^ we obtained a phase-connected solution for 
the pulsar accounting for every turn of the neutron star. In 
short order it became apparent that this pulsar has a large pe- 
riod derivative and therefore a very small characteristic age, 
Tc = P/2P = 7300 yr Some 18 months after we began timing 
the pulsar, on about MJD 52000, it underwent a very large 
rotational glitch, with a fractional period spin-up of two parts 
in one million (see Fig. ^ for the evolution of spin parame- 
ters over time). In magnitude, this is typical of the glitches 
that the « lOkyr-old Vela pulsar experiences every few years 
but, interestingly, there is no evidence for an exponentially 
decaying component as is observed in the Vela glitches (e.g., 
Dodson et al. 2002). 

It is not possible to track the rotation phase of the neu- 
tron star across this large glitch. We find that fitting a tim- 
ing model to the full three-year data set since the glitch shows 
a large quasi-periodicity in the residuals with a ^ 500 day 

* See|http://www.atnf.csiro.au/research/pulsar/tempo|. 
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period, probably due to "timing noise". This can be largely 
absorbed with a fit requiring a total of six frequency deriva- 
tives to whiten the data. For simplicity and to aid observers, 
the solution we present in Tabled is based on the most re- 
cent one year's worth of data, requiring only two derivatives, 
with the second frequency-derivative representing the amount 
of timing noise. We also provide the net jump in spin param- 
eters at the glitch epoch. We emphasize that the value of v is 
not deterministic (see also Fig.^. It yields an apparent brak- 
ing index of rotation n = vv/iv)^ « 40 (where n is defined by 
V ex -v"), which is about 15 times larger than the expected 
long-term deterministic value. 

The presence of significant amounts of timing noise in 
PSR J 1357-6429 is a presumed consequence of the less-than- 
perfect rotational stability of the young, hot neutron star with 
a complex superfluid interior (e.g., Bildsten & Ecstein 1983; 
JArzoumanian et al. 1994; D'Alessandro & McCulloch 1997). 
Such noise biases the pulsar position obtained through tim- 
ing methods. However, we did not have to rely on tim- 
ing measurements to obtain the position, since we performed 
a pulse-gated interferometric observation with the Australia 
Telescope Compact AiTay (ATCA) on 2000 August 29. This 
5-hour observation used the array in its 6A configuration, giv- 
ing a maximum baseline of 6km, with a central frequency of 
1376 MHz and a bandwidth of 128 MHz. The interferometric 
pulsar position, determined from an image made by subtrac- 
tion of an off -pulse data set from an on-pulse data set, is given 
in Tabled and was used for all timing fits. 

3. THE AGE, DISTANCE, AND VICINITY OF PSR J 1357-6429 

The actual age of PSR J 1357-6429, assuming a constant 
magnetic moment, is t = 2Tc{\-{v/vo)"~^'\/{n- 1), where V() 
is the initial rotation frequency of the neutron star. We do not 
know the actual braking index of the pulsar However, if it lies 
in the range 2 < « < 3, as is the case for all four pulsars with n 
measured via phase-coherent timing (see Camilo et al. 200fl 
and references therein), then the real age is still r < 15 kyr and 
may be < 10 kyr, depending on the actual values of n and vq. 
Also, while such youth is not require d by the occurrence of the 
large glitch (e.g.. lHobbs et al.l l2002). it is certainly consistent 
with it. There is therefore no doubt that PSR J 1357-6429 is 
a very young pulsar, one of the ^10 youngest known in the 
Galaxy (see Table|3. 

PSR J 1357-6429 is also apparently nearby: according to 
the ICordes & Lazio' (2002') model for the Galactic free elec- 
tron distribution and our dispersion m easure, the distance 
is t/ = 2.5 kpc. The iTavlor & CordesI (1993) model yields 
d ^ A kpc, suggesting the order of the uncertainty inherent in 
such estimates. The implied distance of the pulsar away from 
the Galactic plane is 100-200pc. In any case, this pulsar ap- 
pears to be one of the two or three nearest to the Earth among 
those known with r < 10 kyr. Considering, furthermore, that 
all known pulsars apparently younger than PSR J 1357-6429 
are associated with a pulsar wind nebula (PWN) or a super- 
nova remnant (SNR), it seems reasonable to expect that this 
might also be the case with PSR J 1357-6429. 

Figure|2]shows the area around PSR J 1357-6429 as seen in 
the 2.4GHz continuum survey of Duncan et al. ( 1995), with 
two SNR candidates proposed by Duncan et al. ( 1997) indi- 
cated. It is not clear what, if any, is the relationship between 
both SNR candidates, or of either one with a much larger su- 
perimposed shell SNR candidate, G3 10.5-3.5. Also, no spec- 
tral, distance or age information is available for these candi- 
dates. The G309.8-2.6 SNR candidate appears at this resolu- 



tion to be about 30' in length with a relatively bright resolved 
"core" to the southwest. In turn, PSR J 1357-6429 is located, 
at least in projection, slightly to the east of this core, but co- 
incident with emission from the SNR candidate. Apart from 
PSR J 1 357-6429, there are seven other known pulsars located 
in the area of the Figure, all old, with 0.5 < Tc < 70 Myr. 

4. DISCUSSION 

PSR J 1357-6429, with r^ = 7300 yr, is among the youngest 
pulsars known in the Galaxy. It is also the second youngest 
pulsar discovered in the Parkes multibeam survey, after 
PSR Jl 119-6127 with r^ = 1600 yr (Camilo et al. 2000). The 
pulsar is located near the SNR candidate G309. 8-2.6. While 
all apparently younger pulsars are associated either with a 
PWN or SNR, three slightly older pulsars are not obviously 
associated with either such object (Table|2j. Further work is 
therefore needed in order to investigate the nature of G309.8- 
2.6, and also to determine whether it is in fact associated with 
the pulsar. 

Being a young and energetic pulsar located at apparently 
only ~ 3 kpc, one would expect substantial X-ray emission to 
be detectable from PSR J 1357-6429. As shown in Table 13 
this is the case for all younger pulsars (interestingly, only the 
Crab pulsar among these is a known optical emitter, and only 
the Crab and PSR B 1509-58 are known 7-ray emitters). A 
check of the HEASARC archives reveals no useful data to 
address this question, and X-ray observations of this source 
will have to be carried out. Also, there is no known EGRET 7- 
ray so urce coincident with this relatively energetic and nearby 
pulsar JHartman e t al. 1999), but PSR J1357-6429 is likely to 
be a good target for the future GLAST mission. 

The large value of i' (and hence apparent n) implied by 
the increasing i) (Fig. ^ possibly results from decay af- 
ter previous (uns een) glitches ( Johnston & Gallowav 199^ 
IWang et al.'2000). However, it is notable that there is no dis- 
cernible jump in !>■ at the time of the observed glitch. In most 
cases, following a large glitch, a portion of the jump in i' and 
i) decays with a time scale that ranges from a few days (for the 
Crab pulsar) to many months (in the Vela pulsar). This results 
in a jump in j>' at the time of the glitch. The absence of this 
for the present glitch suggests either that the decay time-scale 
is much longer than our data span, or that different glitches in 
this pulsar have different characteristics. 

The radio luminosity of PSR J 1357-6429 is L1400 = 
51400^^^ ^ 3mJykpc^, where Shoo is the flux density in the 
1400 MHz band. This is a rather low value: for the 10 radio 
pulsars in Table|2j the median L^oo ^ 25 mJy kpc^. However, 
it is in keeping with the realization that the very youngest 
pulsars (t < 10 kyr) are not obviously more luminous than 
slightly older pulsars (the median luminosity for pulsars with 
Tc < 10^ yr is L1400 ^ 30mJykpc^; Camilo 2004), or than 
middle-aged pulsars (e.g.. Camilo et al. 2002b). 

As we now show, the discovery of a low-luminosity pulsar 
such as PSR J 1357-6429 contributes greatly toward the esti- 
mated birth rate of very young pulsars. To do this, we make 
use of the pulsars listed in Table |2] that were detected in the 
Parkes multibeam survey: B1509-58, Jl 119-6127, J1357- 
6429, B1610-50 and J1734-3333. For each pulsar we carried 
out a Monte Carlo simulation to correct for the volume of the 
Galaxy probed by the multibeam survey, using the V/Vmax- 
style technique to estimate the total number of similar ob- 
jects b eaming toward Earth (for details see, e.g..Lorim eret alJ 
II993L and references therein). Using the most recent esti- 
mates of the Galactic pulsar distribution in these simulations 
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jLorimeil2004h . we find that the total number of such pulsars 
with Tc < lOkyr is 107 ±95. Taking into account the correc- 
tion due to the unknown "beaming fraction" /, the implied 
birth rate of this populati on is S = ( 1.1 ± 1.0) //century"' 
(where / « 0.5; e.g., Frail & Moffett 1993). While having 
a large uncertainty (reflecting the inclusion of PSR J 1357- 
6429, which accounts for 90% of the V/Vmax correction so 
that it dominates the calculation; more generally, we may not 
know enough about the low end of the radio luminosity dis- 
tribution), this birth rate is comparable to those obtained us- 
ing much larger (and on average much older) samples (e.g., 
iLprimer et al. 1993; Lyne et al. 1998). 

Despite the poorly known pulsar birth rate and beaming 
fraction, it appears unlikely that more than « 100 Galactic 
pulsars with t^ < lOkyr are beaming toward us, since the esti- 
mated core-collapse SN rate in our Galaxy (while also having 
significant uncertainties) must account als o for other branches 
of neutron star production (see, e.g., iBrazier & Johnston! 
[B)98; Gaensler et al. 1999, 2000, for discussions of birfli rates 
of other types of neutron stars). In this case, the 1 1 pulsars 
listed in Table |2]located on the "near side" of the Galaxy (all 
but PSR J 1846-0258) could represent about one quarter of 
all such young pulsars from which we may ever detect pulsa- 
tions. Alternatively, the discovery of many more would have 
significant implications for some of the assumptions underly- 
ing these estimates. Nine of the 12 pulsars listed in Table |2l 



have been discovered in the past seven years, in both directed 
and undirected searches, and at radio wavelengths as well as 
in X-rays, methods that suffer from significantly different se- 
lection effects. In view of the above discussion, it is important 
to continue with such searches, employing diverse methods. 

Young pulsars are rare, and young nearby ones are rarer 
still. If PSR J1357-6429 is indeed located at <i < 3 kpc, with 
a visual extinction Ay < 6, and with a probable real age of 
5-15 kyr, it is well possible that its birth event provided a 
spectacular sight to some of our recent pre-historical ances- 
tors. Of greater relevance for us, if not quite as spectacularly, 
future study of this very young and nearby pulsar and any 
possible PWN/SNR companion may contribute one more sig- 
nificant piece toward understanding young neutron stars and 
their environments. 
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TABLE 1 
Measured AND derived parameters for PSR J1357-6429. 



Parameter 

Right ascension (J2000) 

Declination (J2000) 

Rotation frequency, v (s"' ) 

Frequency derivative, v (s^^) 

Second frequency derivative, ;>' (s^^ ) 

Spin period, P (s) 

Period derivative, P 

Epoch (MJD) 

R.M.S. residual (ms) (whitened) 

Epoch of glitch (MJD) 

Frequency step at glitch, t^v (s"') 

Change in v at glitch. A;/ (s^^) 

Dispersion measure, DM (cm^^ pc) 

Flux density at 1374 MHz, S1400 (mJy) 

Pulse profile FWHM, wjo (ms) 

Pulse width at 10% of maximum amplitude, vi'io (ms) 

Surface magnetic field, B (gauss) 

Chai'acteristic age, Tc (kyr) 

Spin-down luminosity, E (ergs s"' ) 

Distance, d (kpc) 

Radio luminosity at 1400 MHz, S\\mS- (mjykpc^). . 



Value 



13'i57"'02?43(2)'' 

-64° 29' 30" 2(1)" 

6.0201677726(4) 

-1.305395(4) X lO"" 

1.16(2) X 10-21 b 

0.16610832750(1) 

3.60184(1) X IQ-'^ 

52921.0 

0.8 

52021 ± 16 

1.46(1) X 10-5 c 

-7.02(3) X lO-"* " 

127.2(5) 

0.44(5) 

15 

31 



7.8 X 10'^ 

7.3 

3.1 X lO^i^ 

~2.5 
~2.7 



Note. — Figures in parentheses (except those for position) are twice the nominal 1 a TEMPO uncertainties in the least-significant digits quoted. The timing 
fit is based on data collected over MJD 52737-53104. 



"Position is obtained from interferometric ATCA observations. 

''This parameter is not stationary (see §|2|for details; see also Fig. ITT. 

■^Assumed glitch epoch for this fit is MJD 52021. The error in Ai^ is dominated by the uncertainty in this epoch. 



TABLE 2 
Parameters for all Galactic pulsars known with apparent age < 10 kyr. 



Pulsar 


SNR/PWN 


P 


log Age" 


logEb 


d 


Wave- 


1 d 

i^.Y pwn 


1 tl 

^.v p.sr 


^1400 




Refs. 






(ms) 


(yr) 


(ergss-') 


(kpc) 


band= 


(10-3£) 


(10-3£) 


(mjykpc^) 






J1846-0258 


Kes75 


323 


2.9 


36.9 


21 


X 


175 


8 






(1) 


J0205-H6449 


3C58 


65 


2.9 


37.4 


3.2 


(x)r 


0.4 


0.007 


~0.5 




(2) 


B0531-H21 


Crab 


33 


3.0 


38.7 


2 


rx 


29 


1 


56 




(3) 


B 1509-58 


MSH 15-52 


150 


3.2 


37.3 


5 


(x)r 


8 


1.7 


35 




(4) 


Jl 119-6127 


G292.2-0.5 


407 


3.2 


36.4 


r^ 6 


rx 


0.02 


0.1 


~29 




(5) 


J1811-1925 


01 1.2-0.3 


64 


3.2 


36.8 


5 


X 


5 


0.6 






(6) 


Jl 124-5916 


G292.0-H1.8 


135 


3.5 


37.1 


>6 


(r)x 


0.2 


0.02 


>3 




(7) 


J1930-H1852 


G54. 1-1-0.3 


136 


3.5 


37.1 


~5 


(r)x 


1.7 


0.3 


~2 




(8) 


J1357-6429 


G309.8-2.6? 


166 


3.9 


36.5 


~2.5 


r 


7 


7 


~3 




(9) 


B1610-50 


— 


231 


3.9 


36.2 


~7.8 


r 


<0.7 


<0.7 


~ 150 




(10) 


J 16 17-5055 





69 


3.9 


37.2 


~6.7 


xr 





1 


~22 




(11) 


J1734-3333 


— 


1169 


3.9 


34.7 


~7.4 


r 


7 


7 


~27 




(12) 


References. - 


— In each given reference, furth 

(3) WeisskoDfetal. i200Q); (4) 


3r references can be found to the X 


-ray observations and dist 


ance estimates, 
Roberts et alj 1 


where appropriate. 


(1) 


Helfande.alJl200.3l,:f2^ 


ISlane et all 12004)1 


Gaensler et al. 


2002); (5) Gonzalez & Safi-Harh 


'2003); (6) 


20031); r7)IHu2hes 


etal 


120031); (S^ICamiloetalJ 



I2QQ2J) : (9) this work; (10) Pivovaroff et al. ■ 2000); (11) Kasoi et al. i 1998); (12) Moms et al. ■ 2002 ). 

Note. — Anomalous X-ray Pulsars and Soft Gamma-ray Repeaters are not included in this list. For definiteness we also do not include a few pulsars that may well have r < 10 kyr, 
but for which r,- > lOkyr, such as the Vela pulsar and PSR J2229-I-61 14 i Halpern et al. 2001). 



^For B0531-I-21, J0205-I-6449, and J1811-1925, historical ages are given assuming associations with SN1054, SN1181, and SN386, respectively IStephenson & Greenl2'Q03) : for 
other pulsars we list the characteristic age, Tc = P/2P. 

Spin-down luminosity, E = A-k IPjP' , where P is the pulsar period and 7=10 g cm". 
'^Pulsar detected in X-rays (x) or radio (r). Discoveiy waveband is listed first, and directed searches are denoted by parentheses. 

X-ray luminosities are all given in the 2-10 keV band and assume isotropic emission. 



THE VERY YOUNG PULSAR J1357-6429 
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Fig. 1 . — The spin history of PSR J1357-6429. The frequency (v) evolution is shown in panels (a) and (b). the latter with an expanded scale after removal 
of u and u fitted to the 100 days just prior to the glitch that occuiTed between MJD 52005 and MJD 52037. The error bars are too small to see in these panels. 
Panel (c) shows the run of i>, obtained from fits at 50-day intervals using approximately 100-day sections of data. The discontinuity in u at the glitch (0.5% of i>) 
reverses about two years' worth of change in u prior to or after the glitch. However, these "continuous" changes in v (the slopes in panel [c]) are not stationary 
but are rather likely due to a stochastic process involving the irregular transfer of angular momentum from the interior to the crust of the neutron star (see §121 
§13 and also, e.g.. .Alpar et al-1984;.Andersson et al-2003.V 
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Fig. 2. — 2.4-GHz image of the Galactic plane in the vicinity of PSR J 1357-6429, using data from the Parkes continuum survey of Duncan et al.'fl995*), at a 
spatial resolution of 10'. The greyscale is linear, and ranges between -0.15 and +1.3 Jy beam"'. The two SNR candidates in this region proposed bv Duncan et all 
119971) are indicated, while the position of the pulsar is marked by a cross (the size of the cross is much larger than the positional uncertainty of the pulsar). 



